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production in polarized proton-proton collisions at y/s = 500 GeV 
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We report the first measurement of the parity violating single-spin asymmetries for midrapid- 
ity decay positrons and electrons from W + and W~ boson production in longitudinally polarized 
proton-proton collisions at ^/s = 500 GeV by the STAR experiment at RHIC. The measured asym- 
metries, Af + = -0.27 ±0.10 (stat.) ±0.02 (syst.) ± 0.03 (norm.) and AY~ = 0.14 ±0.19 (stat.) ± 
0.02 (syst.) ± 0.01 (norm.), are consistent with theory predictions, which are large and of oppo- 
site sign. These predictions are based on polarized quark and antiquark distribution functions 
constrained by polarized DIS measurements. 
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PACS numbers: 14.20.Dh, 13.88.+e, 13.38. Be. 13.85.Qk 



Understanding the spin structure of the nucleon re- 
mains a fundamental challenge in Quantum Chromody- 
namics (QCD). Experimentally, polarized deep- inelastic 
scattering (pDIS) measurements have shown that the 
quark spins account for only w 33% of the proton spin pQ . 
Semi-inclusive pDIS measurements [SHI] are sensitive to 
the quark and antiquark spin contributions separated by 
flavor [5j |6] . They rely on a quantitative understanding 
of the fragmentation of quarks and antiquarks into ob- 
servable final-state hadrons. While the sum of the con- 
tributions from quark and antiquark parton distribution 
functions (PDFs) of the same flavor is well constrained, 
the uncertainties in the polarized antiquark PDFs sepa- 
rated by flavor remain relatively large [6] . 

High-energy polarized proton collisions at y/s = 
200 GeV and y/s = 500 GeV at RHIC provide a unique 
way to probe the proton spin structure and dynamics 
using hard scattering processes [7]. The production of 
W bosons at y/s — 500 GeV provides an ideal tool to 
study the spin-flavor structure of sea quarks inside the 
proton. W^' bosons are dominantly produced through 
U + d (d + u) interactions and can be detected through 
their leptonic decays [8]. Quark and antiquark polarized 
PDFs are probed directly in calculable leptonic W de- 
cays at large scales set by the mass of the W boson. 
The production of W bosons in polarized proton col- 
lisions allows for the observation of purely weak inter- 
actions, giving rise to large parity-violating longitudinal 
single-spin asymmetries. A theoretical framework has 
been developed to describe inclusive lepton production, 
p + p — > W + X —> P + X, that can be directly com- 
pared with experimental measurements using constraints 
on the transverse energy and pseudorapidity of the final- 
state leptons PHO]. 

In this letter, we report the first measurement of the 
parity violating single-spin asymmetries for midrapidity 
decay positrons and electrons from W + and W~ boson 
production in longitudinally polarized p + p collisions at 
y/s = 500 GeV by the STAR experiment at RHIC. The 
asymmetry is defined as Al = (c + — o~)/{o + + a~), 
where is the cross section when the helicity of the 

polarized proton beam is positive (negative). 

The STAR detector systems [TT] used in this mea- 
surement are the Time Projection Chamber [T2] (TPC) 
and the Barrel [T3] and Endcap [T3] Electromagnetic 
Calorimeters (BEMC, EEMC). The TPC provides 
tracking for charged particles in a 0.5 T solenoidal 
magnetic field for pseudorapidities < 1.3 with 
full azimuthal coverage. The BEMC and EEMC are 
lead-scintillator sampling calorimeters providing full az- 
imuthal coverage for |ry| < 1 and 1.09 < rj < 2, respec- 
tively. 

The data analyzed in this letter were collected in 2009 




FIG. 1: (a) Ratios of with respect to the 4x4 BEMC 
E T sum, E^/E^ =4x4 , (b) the near-cone BEMC, EEMC 
and TPC E T sum, E%/ E* =R<0 - 7 , and (c) correlation of 
the signed pr balance variable and E^p. MC shape distribu- 
tions (arbitrary normalization) are shown in (a) and (b) for 
W ± — > e + X as filled histograms in comparison to both data 
distributions. 



with colliding polarized proton beams at s/s = 500 GeV 
and an average luminosity of 55 x 10 30 cm _2 s _1 . The 
polarization of each beam was measured using Coulomb- 
Nuclear Interference proton-carbon polarimeters [15] . 
which were calibrated using a polarized hydrogen gas-jet 
target [16 . Longitudinal polarization of proton beams 
in the STAR interaction region (IR) was achieved by 
spin rotator magnets upstream and downstream of the 
IR that changed the proton spin orientation from its sta- 
ble vertical direction to longitudinal. Non-longitudinal 
beam polarization components were continuously moni- 
tored with a local polarimeter system at STAR based on 
the Zero-Degree Calorimeters with an upper limit on the 
relative contribution of 15% for both polarized proton 
beams. The longitudinal beam polarizations averaged 
over all runs were Pi = 0.38 and P2 = 0.40 with cor- 
related relative uncertainties of 8.3% and 12.1%, respec- 
tively. Their sum Pi + P2 — 0.78 is used in the analysis 
and has a relative uncertainty of 9.2%. 

Positrons (e + ) and electrons (e~) from W + and W~~ 
boson production with |?y e | < 1 are selected for this anal- 
ysis. High-pT e 1 *" are charge-separated using the STAR 
TPC. The BEMC is used to measure the transverse en- 
ergy E?p of e + and e~ . The suppression of the QCD back- 
ground is achieved with the TPC, BEMC, and EEMC. 

The selection of W candidate events is based on kine- 
matic and topological differences between leptonic 
decays and QCD background events. Events from 
decays contain a nearly isolated e ± with a neutrino in 
the opposite direction in azimuth. The neutrino escapes 
detection leading to a large missing energy. Such events 
exhibit a large imbalance in the vector px sum of all re- 
constructed final-state objects. In contrast, QCD events, 
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Lepton charge sign x l/p T (c/GeV) 

FIG. 2: Et as a function of the ratio of the TPC recon- 
structed charge sign to the transverse momentum pr- The 
black-solid and red-dashed lines indicate the selected kinematic 
region used for the asymmetry analysis. 

e.g. dijet events, are characterized by a small magnitude 
of this vector sum imbalance. 

Candidate W events were selected online by a two-step 
energy requirement in the BEMC. Electrons or positrons 
from W production at midrapidity are characterized by 
large E T peaked at « M w /2 (Jacobian peak). At the 
hardware trigger level, a high tower calorimetric trig- 
ger condition required Et > 7.3 GeV in a single BEMC 
tower. At the software trigger level, a dedicated trig- 
ger algorithm searched for a seed tower of Et > 5 GeV 
and computed all four possible combinations of the 2x2 
tower cluster E T sums and required at least one to be 
above 13 GeV. A total of 1.4 x 10 6 events were recorded 
for a data sample of 12pb _1 . A Vernier scan was used 
to determine the absolute luminosity [17] . 

An electron or positron candidate is defined to be any 
TPC track with p T > lOGeV/c that is associated with 
a primary vertex with \z\ < 100 cm, where z is measured 
along the beam direction. A 2 x 2 BEMC tower cluster 
Et sum E T , whose centroid is within 7 cm of the pro- 
jected TPC track, is required to be larger than 15 GeV. 
The excess BEMC Et sum in a 4 x 4 tower cluster cen- 
tered around the 2x2 tower cluster is required to be 
below 5%, as indicated by the vertical dashed line in 
Fig. [l] (a). A cone, referred to as the near-side cone, is 
formed around the e ± candidate with a radius R = 0.7 
in r}-(j> space. The excess BEMC, EEMC, and TPC E T 
sum in this cone is required to be less than 12% of the 
2x2 cluster Et, as shown in Fig. [I] (b) by the vertical 
dashed line, pythia 6.205 [IS] Monte-Carlo (MC) shape 
distributions (arbitrary normalization) for W ± — > e + v 
passed through the geant [19] model of the STAR detec- 
tor are shown in Figs. [I] (a) and[T](b) as filled histograms 
motivating both ratio cuts. The missing energy require- 
ment is enforced by a cut on the pt balance vector, de- 
fined as the vector sum of the e candidate pt and the 
Pt vectors of all reconstructed jets, where the jet thrust 
axis is required to be outside the near-side cone. Jets 
are reconstructed using a standard mid-point cone algo- 
rithm used in STAR jet measurements J2U] based on the 



TPC, BEMC, and EEMC. A scalar signed p T balance 
variable is formed, given by the magnitude of the pt bal- 
ance vector and the sign of the dot-product of the pt 
balance vector and the electron pt vector. This quantity 
is required to be larger than 15 GeV/c. The correlation 
of the signed pt balance variable and E T is shown in 
Fig. [I] (c) . The range for accepted W candidate events is 
marked by red dashed lines. The lower cut is chosen to 
suppress the contribution of background events whereas 
the upper cut is mainly applied to ensure proper charge 
sign reconstruction. Background events from Z° — > e + e~ 
decays are suppressed by rejecting events with an addi- 
tional electron/positron-like 2x2 cluster in the recon- 
structed jet where the E T x2 > p J ^ t /2 and the invariant 
mass of the two electron/positron-like clusters is within 
70 to 140GeV/c 2 . This avoids Z° contamination in the 
data-driven QCD background described below. 

Figure [2] shows E T as a function of the ratio of the TPC 
reconstructed charge sign to the transverse momentum 
Pt for electron and positron candidates that pass all the 
cuts described above. Two well-separated regions for pos- 
itive (negative) charges are visible, identifying the W +Iy ~' 1 
candidate events up to E T « 50 GeV. The range of E T 
for accepted W candidate events, 25 < E T < 50 GeV, 
is marked by red dashed lines. Entries outside the black 
solid lines in Fig. [2] were rejected due to false track re- 
construction. 

Figure [3] presents the charge-separated lepton E T dis- 
tributions based on the selection criteria given above. 
W candidate events are shown as the black histograms, 
where the characteristic Jacobian peak can be seen at 
~ Mw/2. The total number of candidate events for 
is 462(139) for 25 < E T < 50 GeV indicated 
by vertical red dashed lines in Fig. [3] The number of 
background events was estimated through a combination 
of pythia 6.205 [H] MC simulations and a data-driven 
procedure. The e + ^~' > background from W +t "~^ boson 
induced t + (' decays and Z° — > e+ + er decays was es- 
timated using MC simulations tobel0.4±2.8(0.7±0.7) 
events and 8.5 ±2.0 events (identical for both e + ( - )), re- 
spectively. The remaining background is mostly due to 
QCD dijet events where one of the jets missed the STAR 
acceptance. We have developed a data-driven procedure 
to evaluate this type of background. We excluded the 
EEMC (1.09 < r) < 2) as an active detector in our analy- 
sis to estimate the background due to missing calorimeter 
coverage for —2 < rj < —1.09. The background contri- 
bution due to missing calorimeter coverage along with r 
and Z° background contributions have been subtracted 
from both W +(i -~^ E T distributions. The remaining back- 
ground, presumably due to missing jets outside the STAR 
1 77 1 < 2 window, is evaluated based on an extrapolation 
from the region of E T < 19 GeV in both E T dis- 

tributions. The shape is determined from the E T distri- 
bution in events previously rejected as background with 
systematic variations of the signed pt balance cut be- 
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FIG. 3: Et for W + (bottom) and W~ (top) events showing 
the candidate histograms in black, the full background esti- 
mates in blue and the signal distributions in yellow. 



low 15GeV/c. This shape Ej, distribution is normalized 
to both E%, distributions for Ef, < 19GeV. The 

total number of background events for e + '~) is 39 ± 9 
(23 ± 6) for 25 < E% < 50 GeV shown in Fig. [| as the 
blue line. The errors on the total background are mostly 
from the data-driven background events. 

The leptonic asymmetry from W ± decay, AY , was 
obtained from: 



aY ± = 



1 2 gj 
3± Pi + P 2 ~ 



Y> i R,N} v± 



(1) 



where P\^, are the mean polarizations, are 
candidate yields for all four beam helicity configura- 
tions i = {++,H — , — hj }, and Ri are the respec- 
tive relative luminosities. The longitudinal single-spin 
asymmetry At, for Z° bosons has been estimated us- 
ing a full next-to- leading (NLO) order framework [TP] . 
With the W ± selection criteria we estimated the Z° 
asymmetry to be Af = —0.06. This value has been 
used to determine the polarized background contribution 
a+C - ) = -0.002 ± 0.001 (-0.005 ± 0.002). The unpo- 
larized background correction for W ± candidate events 
is /?+(-) = 0.938 ± 0.017(0.838 ± 0.032). This dilution 
factor is due to background events passing all W selec- 
tion cuts and is determined by (3 = S/(S + B), where 
S (B) is the number of signal (background) events for 
25 < E%, < 50 GeV. 

The relative luminosities Ri — ^ fc Mfc/(4Mj) are de- 
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termined from the ratios of yields Mj of QCD events, for 
which parity conservation is expected. The M, are statis- 
tically independent from because the isolation cut 
on the 2x2 / 4 x 4 tower Ex sum, shown in Fig. [I] was 
reversed for those events. Additionally, an upper limit of 
20 GeV was set on E T . 

Figure [4] shows the measured leptonic asymmetries 
-0.27 ± 0.10 (stat.) ± 0.02 (syst.) and = 
0.14 ± 0.19 (stat.) ± 0.02 (syst.) for \r/ e \ < 1 and 25 < 
Ej. < 50 GeV. The vertical black error bars include only 
the statistical uncertainties. The systematic uncertain- 
ties are indicated as grey bands. The statistical uncer- 
tainties dominate over the systematic uncertainties. The 
asymmetry A^ observed in statistically independent sam- 
ples of QCD dominated events was found to be 0.04±0.03 
(0.00 ±0.04) for positive (negative) charged tracks and is 
consistent with zero. We assumed the experimental limit 
on the polarized background Al to be 0.02 as a system- 
atic uncertainty of AY ■ This limit on polarized back- 
ground and the uncertainty in unpolarized background 
dilution have been added in quadrature to account for 
the total systematic uncertainty of AY ■ The normal- 
ization uncertainty of the measured asymmetries due to 
the uncertainty for the polarization sum Pj + P 2 is 0.03 
(0.01) for AY +{ ■ The normalization uncertainty is of 
similar size as the systematic uncertainty of the asymme- 
try measurement. 

In Fig. [4j the measured asymmetries are compared to 
predictions based on full resummed (rhicbos) [5] and 
NLO (che) [10] calculations. The che calculations use 
the DSSV08 polarized PDFs [5], whereas the RHICBOS 
calculations are shown in addition for the older DNS-K 
and DNS-KKP [3T] PDFs. The che and rhicbos re- 
sults are in good agreement. The range spanned by the 
DNS-K and DNS-KKP distributions for Ad and Au co- 
incides, approximately, with the corresponding DSSV08 
uncertainty estimates [51 [B]. The spread of predictions 
for AY +i ' is largest at forward (backward) r\ e and is 
strongly correlated to the one found for the d (u) polar- 
ized PDFs in the RHIC kinematic region in contrast to 
the backward (forward) r/ e region dominated by the be- 
havior of the well-known valence u (d) polarized PDFs 
|10) . At midrapidity, production probes a combi- 

nation of the polarization of the u and d (d and u) quarks, 
and AY +i ' is expected to be negative (positive) [5j[6]. 
The measured AY + is indeed negative stressing the di- 
rect connection to the u quark polarization. The central 
value of AY is positive as expected with a larger sta- 
tistical uncertainty. Our At, results are consistent with 
predictions using polarized quark and antiquark PDFs 
constrained by inclusive and semi-inclusive pDIS mea- 
surements, as expected from the universality of polarized 
PDFs. An independent measurement of W boson pro- 
duction from RHIC is being reported by the PHENIX 
collaboration 122 . 
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In summary, we report the first measurement of the 
parity violating single-spin asymmetries for midrapid- 
ity |j7 e | < I, decay positrons and electrons from W + 
and W~ boson production in longitudinally polarized 
p + p collisions at s/s = 500 GeV by the STAR exper- 
iment at RHIC. This measurement establishes a new 
and direct way to explore the spin structure of the pro- 
ton using parity-violating weak interactions in polarized 
p + p collisions. The measured asymmetries probe the 
polarized PDFs at much larger scales than in previous 
and ongoing pDIS experiments and agree well with NLO 
and resummed calculations using the polarized PDFs of 
DSSV08. Future high-statistics measurements at midra- 
pidity together with measurements at forward and back- 
ward pseudorapidities will focus on constraining the po- 
larization of d and u quarks. 
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